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ABSTRACT 


A  solution  CO  the  problem  of  strong  Interaction  between  the 
shock  wave  and  the  boundary  layer  has  been  obtained  for  the  case  where 
velocity  slip  and  temperature  jump  boundary  conditions  are  consistent 
at  the  wall.  It  Is  shown  that  the  addition  of  slip  boundary  conditions 
yields  a  correction  of  order  \jS/X}  to  the  no  slip  solution. 

Estimates  are  made  of  the  effect  of  slip  on  Induced  pressures 
and  skin  friction  for  the  case  of  the  adiabatic  wall.  In  addition,  It 
Is  shown  that  the  Inclusion  of  slip  boundary  conditions  does  not  change 
the  energy  transfer  to  the  wall  from  the  no  slip  values. 
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STRONG  INTEBACTION  WITH  SLIP  BOUNDARY  CONDITIONS 


J.  Aroesty 


1  INTRODUCTION 

The  status  of  the  hypersonic  leading  edge  problem  has  not  changed 
appreciably  in  the  period  since  the  earliest  recognition  of  the  shock  wave¬ 
boundary  layer  interaction  phenomenon.  The  theoretical  aspects  of  the  problem 

have  mainly  been  considered  from  the  viewpoint  of  boundary  layer  theory,  with 

fl23l  CaI 

few  exceptions^  *  *  .  There  has  been  some  work  by  Charwat^  on  the 

application  of  a  near-free-molecule  technique  to  the  region  very  close  to 

the  tip.  This  approach,  which  actiially  considers  the  various  types  of 

collisions  encountered  by  molecules  which  have  been  reflected  off  the  plate 

surface,  is  similar  to  the  "physical  iteration"  scheme  used  by  Charwat  and 

Baker^^^  in  their  consideration  of  near-free-molecule-f low  over  a  sphere. 

This  counting  of  collisions  becomes  quite  complicated  as  more  and  more 

collisions  occur  and  it  becomes  increasingly  difficult  to  assign  a  mean 

free  path  to  each  type  of  molecular  encounter. 

Another  attack  on  the  problem,  not  within  the  scope  of  boundary 

[6 1 

layer  theory,  was  carried  out  by  Laurmann^  who  considered  the  problem  of 
slip  flow  over  a  short  flat  plate.  The  solution  was  based  on  a  linearization 
of  the  full  Navier-Stokes  equations  about  the  free  stream  velocity.  The 
interesting  possibility  was  raised  that  for  this  problem,  with  sufficient 
slip  at  the  surface,  the  Oseen  solution  could  be  valid  over  the  entire  short 
plate.  Although  it  is  difficult  to  relate  this  problem  of  the  short  plate 
to  that  of  the  leading  edge  region  of  a  long  flat  plate  on  physical  grounds, 
Laurmann' result  of  a  near  linear  displacement  thickness  is  certainly 
suggestive  of  the  possible  effects  of  slip  boundary  conditions  at  the  surface. 
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Except  for  a  recent  paper  by  Ogucht^  the  remainder  of  the  analytical 


attempts  using  the  continuum  approach  have  been  within  the  framework  of  hypersonic 

fg] 

similar  solutions  of  the  boundary  layer  equations.  The  earliest  solutions 
suggested  the  possibility  of  a  strong-interaction  region,  where  the  pressure 
grows  like  x  and  the  boundary  layer  grows  like  x^^^.  The  important 


parameter  in  this  region  was  shown  to  be  H 
given  as  functions  of  this  hypersonic  interaction  parameter,  yt  . 


,  and  all  solutions  were 


Although  the  early  measurements  of  the  pressure  distributions  on  the 

sharp  leading  edge  flat  plate  in  hypersonic  flow  seemed  to  confirm  the  existence 

of  this  strong  interaction  region,  later  work  done  at  Berkeley  and  at  Schenectady 

demonstrated  clearly  that  this  parameter  alone  is  certainly  not  able  to  describe 

the  pressures  as  the  plate  tip  is  approached.  Starting  from  the  downstream  side, 

it  becomes  apparent  that  the  pressures  increase,  at  least  qualitatively,  as 

predicted  by  the  strong  interaction  results  until  they  reach  a  plateau  of 

constant  pressure,  and  then  decrease.  The  level  of  the  pressure  plateau  is  a 

function  of  Mach  number  for  sufficiently  large  values  of  .  Since  the 

various  data  have  been  taken  at  a  variety  of  wall  temperatures,  it  is  simplest 

to  characterize  the  experiments  done  at  Berkeley  as  the  insulated  case,  and 

the  experiments  done  at  General  Electric  Company  as  the  highly  cooled  case, 

Even  the  strong  interaction  theories  suggest  that  the  induced  pressures 

increase  with  increased  T so  that  it  is  not  surprising  that  the  level  of 

this  pressure  plateau  is  also  a  function  of  this  temperature  ratio.  Examination 

of  the  available  data  suggests  the  possibility  that  the  extent  of  the  strong 

_l/2 

interaction  region,  where  the  pressure  grows  like  x  '  ,  is  probably  quite 
small,  if  it  exists  at  all. 

Prior  to  the  publication  of  the  data  which  showed  the  existence  of 
a  pressure  plateau  near  the  tip,  attempts  to  modify  the  early  strong  interaction 


2 


solutions  vere  made  by  Oguchi  and  Stewartson.^^^’ ^^^using  hypersonic  small 
disturbance  theory.  On  purely  analytical  grounds,  these  more  recent  solutions 
are  attractive  in  that  they  provide  a  complete  solution  to  the  flow  field 
between  the  shock  wave  and  the  boundary  layer,  whereas  the  early  work,  which 
utilized  the  tangent  wedge  formula,  gave  no  Information  at  all  conceiming  the 
nature  of  the  flow  field  external  to  the  boundary  layer. 


2  OGUCHI 'S  AHALYSES 

Recently  Oguchi gave  a  more  complete  Joining  of  the  boiindary 
layer  and  the  invlscld  flow,  where  the  vorticlty  and  temperature  at  the  edge 

of  the  boundary  layer  are  no  longer  taken  to  be  zero,  as  In  the  earlier 

« 

solutions,  but  rather  at  the  values  that  obtain  from  the  solution  of  the  zero 
order  Invlscld  equations.  Thus,  the  complete  similarity  of  the  flow  field 
enables  one  to  obtain  Information  which  Is  not  at  all  available  In  the  tangent 
wedge  approach,  which  yields  p/p^  alone  at  the  edge  of  the  boundary  layer. 

However,  the  means  necessary  to  achieve  complete  similarity  is  to 
assume  that  the  shock  wave  grows  like  Ax*^,  where  1  >  n  £  l/3.  With  the 
application  of  the  approximate  conservation  laws  across  the  shock,  this  implies 
that  the  pressure  and  temperature  approach  a  singularity  as  the  leading  edge  is 
approached.  Thus,  this  simple  form  of  the  shock  shape,  essential  to  a  similar 
invlscld  flow  field,  is  then  responsible  for  the  high  pressures  near  the  leading 
edge,  and  also  for  the  perhaps  spuriously  hot  gas  at  the  edge  of  the  boundary 
layer. 


The  matching  between  the  shock  layer  and  the  boundary  layer  problem 


Is  effected  by  an  expansion  in  a  small  parameter. 


,  which 


is  suggested  by  the  forms  of  the  temperature  and  vorticlty  variation  at  the 


3 


edge  of  the  viscous  regions.  For  air,  the  small  parameter  is  effectively  of 
order  ^  ,  which  Introduces  some  question  concemisg  its  validity  in  a 

boundary  layer  solution.  The  additional  consideration  of  the  vorticity  due 
to  the  curved  leading  edge  shock  wave  yields  pressure  levels  which  are  higher 
even  than  those  of  the  earlier  solutions. 

In  the  light  of  experiments  that  seem  to  point  to  a  pressure  plateau 

as  the  leading  edge  is  approached,  it  is  possible  that  the  leading  edge  shock 

wave  is  not  very  highly  curved,  and  in  fact,  it  is  certainly  possible  that  it 

is  nearly  straight.  With  this  in  mind,  Oguchi  its  recently  attempted  to 

settle  this  problem  of  shock-boundary  layer  interaction  in  a  manner  which  is 

similar  to  the  viscous  layer  problem  on  a  blunt  body.  Oguchi  has  assumed 

that  the  flow  in  the  entire  region  between  the  shock  wave  and  the  body  may 

be  represented  by  the  boundary  layer  equations.  This  implies  that  the 

entire  pressure  field  in  the  very  leading  edge  region  is  constant,  since  the 

assumption  is  also  made  that  the  shock  wave  is  straight.  In  addition  to  the 

usual  condition  that  u  =  u^  at  the  shock,  it  is  also  required  chat  the 

mass  flow  entering  the  shock  wave  is  equal  to  the  mass  flow  in  the  boundary 

layer.  This  is  done  by  equating  the  boundary  layer  stream  function  at  the 

still  unknown  shock  wave  height,  y  ,  to  the  stream  function  corresponding 

s 

to  the  free  stream  at  that  height. 

Using  the  nomenclature  of  reference  [  8],  this  new  condition  is 


Thtfcyf;)  -  /oe 


(1) 


(2) 


where 

1^03 


•A  n 


(V\ 
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It  is  Oguchl's  hope  that  the  aolutlon  of  the  boundary  layer  equations, 
with  this  additional  boundary  condition, yield  a  zero  value  for  the  shear  stress 
juat  behind  the  shock  wave.  This  la  critical  to  the  Justification  of  the  use  of 
the  boundary  layer  equations  In  the  entire  field,  since  the  shock  wave  is  assumed 
to  be  a  straight  discontinuity,  satisfying  the  appropriate  lUinkine-Hugoniot 
conditions.  Since  the  shock  thickness  varies  inversely  as  the  strength,  It  Is 
more  likely  that  one  Is  able  to  give  an  a  priori  justification  for  the  neglect 
of  the  shock  structure  in  the  case  of  a  very  strong  shock,  such  as  a  normal 
shock,  than  for  this  case  of  a  comparatively  weak  shock. 

In  attempting  to  meet  the  boundary  conditions  for  this  problem, 

Oguchi  utilises  an  assumed  form  for  the  stream  function  which  is  inconsistent 
with  the  resulting  equations,  except  for  the  region  where  <<  /  ,  corres¬ 
ponding  to  the  region  closest  to  the  tip.  ?or  larger  values  of  ,  a 

numerical  solution  is  proposed  which  is  called  a  "local  similarity  solution". 

By  means  of  these  two  solutions,  one  analytic  and  one  numerical,  the  Induced 
pressures  over  the  entire  leading  edge  region  are  predicted,  and  are  shown  to 
agree  quite  well  with  the  more  recent  experiments. 

However,  if  the  value  for  the  sheer  behind  the  edge  of  the  shock 
is  evaluated  properly,  l.e.. 


5 


the  solution  for  the  Limiting  case  of 


yields  the  result  that 


—  J  throughout  the  boundary  layer.  Similarly,  for  the  case 

'X/ 

corresponding  to  the  "locally  similar"  solutions,  ia  always 

Ti*/*ku 

greater  than  zero,  and  in  fact,  is  always  greater  than  0.3.  Thus,  we  see 
that  the  requirement  of  the  vanishing  of  the  shtiar  behind  the  shock  wave  is 
not  met,  and  an  Inconsistency  develops  at  the  shock  wave.  The  way  out  of  this 
inconsistency  is  certainly  not  clear,  but  it  is  apparent  that  the  solution,  as 
it  stands  now.  Is  certainly  open  to  question. 


3  EFFECTS  OF  SLIP 

Laurmann's  suggestion  that  a  linear  displacement  thickness  is  possible 
for  large  slip  velocities,  and  the  well  known  Incompressible  result,  that  small 
slip  at  the  surface  decreases  the  displacement  thickness,  all  point  to  the 
likelihood  that  Inclusion  of  first  order  slip  boundary  conditions  in  the  strong 
Interaction  problem  will  tend  to  reduce  the  curvature  of  the  boundary  layer 
thickness,  and  consequently  reduce  the  pressure  levels  from  the  no  slip  value 
as  the  leading  edge  is  approached.  Street^^'^^  and  Laurmann^ ^  have  considered 
this  problem,  the  latter  using  an  integral  technique,  and  the  former  using  an 
expansion  procedure  somewhat  analogous  to  the  Oguchi  procedure  for  matching 
the  vorticity. 

The  procedure  which  Street  and  Laurmann  used  is  the  expansion  of 

certain  flow  quantities  in  terms  of  a  small  parameter,  which  is  related  to  the 

f  1 3 1 

local  mean  free  path  at  the  wall.  This  technique  is  quite  well  known^^  ,  and 

is  characterized  as  possessing  solutions  which  are  related  to  the  y-derlvatives 

of  the  zero*"^  order  solutions. 

!  1  o  J 

Cassaccio  has  attempted  to  generalize  this  result  to  the  case  of  a 
similar  compressible  flow,  where  no  effect  on  the  external  flow  Is  considered. 
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The  solutions  obtained  are  characterized  by  this  linear  dependence  on  the 

til 

y-derivatlves  of  the  zero'"  '  order  solution. 

If  the  expansion  is  cast  in  the  Howarth-Dorodnits;  n  variables' 
the  appropriate  expansion  parameter  becomes  evident  from  the  expression  for 
at  the  wall,  i.e.,  slip  boundary  conditions. 


f'-  ufo)  --  ‘i>  c  • 


(0) 


(7) 


The  appropriate  expansion  parameter  becomes 


(S) 


Expanding 


H-  ^ 


(S') 


(i.e.,  perturbing  a  similar  solution),  and  inclusion  in  the  momentum  equation 
yields 


zero 


th 


order  I  £ 


(10) 


(11) 
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The  requirement  for  similarity  in  the  2ero^“  order  solution  was  that 


2.  Ut  f 


JL-  =  b/  (13) 

For  the  assumed  first  order  expansion  to  be  consistent,  it  is  clearly  necessary 
that  ^  must  be  of  the  form^JP^  when  K  Is  a  constant.  Since 

is  independent  of  ^  ,  being  determined  only  by  the  wall  temperature,  this 

is  equivalent  to  requiring  that  vary  like  a  power  of  JP  ,  This  is 

Inconsistent  with  the  requirement  for  the  zeroth  order  solution,  which  is  that 

U,  .rtf’ 

)'  ,_j  r.f  f  <“> 

Only  for  the  special  cases  of  ,  corresponding  to  a  flat  plate.  Mi 

corresponding  to  hypersonic  flow,  or  corresponding  to  low  speed  flow, 

f  17] 

is  the  expansion  procedure  completely  consistent.  Maslen*^  gave  the  complete 
solution  tor  the  case  where  ,  corresponding  to  compressible  slip  flow 

over  a  flat  plate.  Otherwise,  It  is  possible  to  show  that  the  only  solution 
satisfying  this  assumed  foru,  for  f  is  the  one  which  drives  the  term  in  the  brackets 
in  equation  (11)  identlally  to  0,  or  .  However,  this  is 

certainly  a  very  special  case,  since  we  have  only  given  one  solution  to  a 


third  order  differential  equation  for  In  general,  we  can  see  that  this 

assumed  expansion  is  not  satisfactory. 
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4  PRESENT  ANALYSIS 


The  present  analysis  will  consist  of  a  perturbation  about  the  strong 
interaction  solution,  in  order  to  represent  the  effects  of  the  slip  boundary 
conditions. 

The  following  assumptions  will  be  made  in  order  to  obtain  the  familiar 
zero  order  equations  given  in  [  3  ],  page  305. 

1.  The  boundary  layer  equations  are  valid. 

2.  Pr  -  1 


N  -  (constant) 


'"Sb  '  U 

-fi,  J>'>l  t  ■  constant 


5.  Tangent  wedge  formula  is  valid  to  describe  the  interaction 
between  the  displacement  thickness  and  the  induced  pressure. 

The  boundary  layer  equations,  with  first  order  slip,  for  Pr  »  1  may  be 
written  as: 


(15) 


(16) 


(17) 


B.C. 


y~C  I  VsO  \  Q  (1®5  Pi.  “  local  mean  free  path 

^  fo 

(19) 

,  >4~ 


(20) 


4*  1  Derivation  of  the  Zero  &  First  Order  Equations 

The  Howarth  Dorodnytain  transformations  are  applied  in  order  to 
render  the  equations  amenable  to  a  similar  solutions  analysis. 


(21) 

(22) 

3  -  % 

(23) 

y  » 

^  (24) 

The  momentum  equation  becomes 

The  energy  equation  becomes 

The  momentum  equation  can  be  put  in  a  somewhat  different  form  for 
the  case  where  the  external  flow  is  hypersonic--  i.e., 

c/f 

f  '  Wj  Jy 


(27) 


\ 

\4,  a.y 
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If  we  set  **  constant,  in  the  Integral  , 


Thus 

Jl,  * 

4>f 


<29  ) 


assuming  a  perfect  gas. 


or  then 

’^T  ^  ^ 


Since  we  are  assuming  an  isentropic  external  flow, 


”  djL 


(31) 


(32) 


or 


A  r 


f-p  </x  ^  77 


-  "  tif 

^  Hv 


(33) 


for  perfect  gas 


for  / 
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Rewriting 


7f>  '? 


for  perfect  gas 
(34) 


^4  /  ?  y  05) 

The  momentum  equation  may  then  he  written 


(36) 


Th.  .Up  ^ 


as 


or 


9(^)  -i  -) 

V^J/y  5  0  ^ 

tj*0 


(37) 


Temperature 

Jump 


or 


^  ^(5)  “ 

T  W 


(38) 


Velocity 

Slip 
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The  form  of  the  slip  boundary  conditions  suggests  that  the  appropriate 
expansion  parameter  for  a  perturbation  solution  about  the  zero  order  strong 


interaction  solution  would  be 


j;  - 


parameter  in  terms  of  M  and  Re  shows  that 
where  ^  [-Zl] 


,  where 


.  Evaluation  of  this  small 

£- 

.  (See  Appendix  A.) 


The  expansion  is  applied  directly  to  the  transformed  equations  — 


•  hi' *  ^ 


At  the  edge  of  the  boundery  layer*  as  17  »  ,  the  following  terms  are 


expanded 


Hi  »  ^ 


where  Pj^,  Uj^,  h^ 
are  constants 


*  hi, + 

Because  the  external  flow  is  adiabatic,  the  stagnation  enthalpy  remains 


conitant  at  Hb> 


When  these  forma  are  placed  in  the  transformed  equations,  the  following 
sets  of  equations  are  obtained 
Zero^^  order  in  £ 


13 


At 


(48) 


At 


4^-»/ 


(45) 


The  •olutions  of  the  zero  order  equations  are  tabulated  in 

,\ 

Ffrst  Order  In  ^ 


and 


^<?'77  •  -  ^iS> 


Hewrltlng 

A/^»tjij»|  ^  a 


n 


■*■  <^7';  •*■  1^7^ 


Where  B.  C.  are 

f'fy)  -  a  (i% 

f,  V*')  •  o 

f.M  -  0 

-  ji/o) « ^o) 
?/*>  -  d 


(53) 


(54) 


(55) 
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Cl8j 


(50) 


(51) 


(52) 


i 


A. 2  Solution 


The  energy  equation  may  be  Integrated  Immediately 


where  and  are  obtained  by  variation  of  parameters. 


Tho  rondltlon  at  W  *  (3  permits  the  constants  A  and  C 
to  be  evaluated,  since  all  the  other  terms  vanish  there. 

However,  the  condition  at  eo  is  identically  satisfied  by  the 
functions,  irrespective  of  the  value  of  the  constants  B  and  D.  It  is 
then  necessary  to  introduce  some  other  requirement  on  the  solution  in  order 
to  establish  a  physically  sensible  solution.  If  the  requirement  is  made  that 
the  various  thicknesses  which  characterize  the  boundary  layer  possess  finite 
dimensions,  it  is  possible  to  evaluate  the  constants  B  and  D,  (See  Appendix 
G). 
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When  these  constants  are  evsluatedj  the  solution  co  the  first  order 


energy  esjuation  becomes 


The  solution  to  the  first  order  momentum  equation  can  be  written 


Since  the  functions  f^,  f^"',  g^'  and  have  already  been  calculated  ^  ^ 

it  would  be  possible  to  evaluate  the  velocity  and  energy  perturbations  as 
numerical  functlms  of  O  . 


In  order  to  calculate  Pj^,  It  is  necessary  to  first  calculate  the 
displacement  thickness,  and  then  utilize  the  tangent  wedge  formula. 

If  we  neglect  terms  of  aftJ  ,  the  dlsplacefflent  thlrkneas  can 


i.'H  J 


be  written 


(62) 


Expanding  in  the  small  parameter  ,  using  the  perfect  gas  law,  and  the 

results  that  Mi"*©  ,  we  may  write 


(63) 
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where 


X  “  /  (64) 


The  Integral,  ,  will  be  in  the  form  ,  since  the  perturbation 

velocity  profile  Is  given  as  a  linear  function  of  p^.  Since  ^ 
may  be  written 

The  use  of  the  tangent  wedge  formula  requires  that  the  slope  of  the 


displacement  thickness  be  evaluated 


if ^  f  7  + 


L 

X 


JA  a  isl 


(67) 


(68) 


(66) 


or 


The  tangent  wedge  formula,  as  given  in  reference  [  8  ]  is 
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Since 


then  b  'A 

^  ^17 


(73) 


The  calculation  of  for  the  general  case  is  quite  complicated, 

since  the  numerical  evaiiiation  of  several  integrals  is  required  in  order  to 
evaluate  f^'  and  hence,  .  Although  the  complete  calculation  is 

possible  in  principle,  it  is  of  questionable  utility,  since  the  expansion 
parameter  is  of  order  .  However,  in  order  to  give  an  estimate  of 


the  magnitude  of  p^,  we  shall  consider  the  case  of  the  adiabatic  wall.  In 

addition,  we  shall  set  P,  ^  ^  ^ 

1=0  whenever  it  appears  in  the  velocity  profile. 

This  is  equivalent  to  setting  the  non-homogeneous  terms  to  zero  in  the 

differential  equations.  If  we  suggest,  on  intuitive  grounds,  that  the  primary 

effect  of  slip  on  the  velocity  profiles  is  not  through  the  change  in  the 

external  flow  field  (i.e.,  p)  but  rather  directly  through  the  change  in  the 

boundary  conditions,  we  can  write  the  approximate  solution  of  the  momentum 

equation  as 


fX)  a 


(74) 


Evaluating  ^  ,  since  g^  ■»  0,  for  the  adiabatic  wall 

f<  *  -  (75) 

or 

*  -ad 

S,  =  -o- 


(76) 


and 


(77) 


‘  ^  So 


(78) 
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From  Oguchi's  paper^  So  ia  evaluated  from  a  solution  of  the  zero*"^  order 

momentum  equation,  Oguchi's  notation  uses  I  to  represent  our  , 


=  1.310  for  air  (79) 


Therefore,  the  nondlmens tonal  coefficient  for  the  pressure  perturbatio-.  is 

Pji  "  -  a.  For  the  usual  case  of  completely  diffuse  surface  interaction,  this 

2  *  6 

constant  a  may  be  replaced  by  1,  since  a  =>  — -7 —  x  2  x  0.499. 


As  demonstrated  in  Appendix  A, 


Therefore,  our  estimate  of  the  induced  pressure  for  tM  adiabatic  wall  is 


(80) 


or 


inr^ 


(81) 


This  can  be  written 


or 


t1 

Sine.  q- 

h 

■  -^3jL  ] 


(82) 


(84) 


An  examination  of  equation  (61)  suggests  thcL  ibis  is  a  i 

th 

the  decrease  in  induced  pressure  from  the  zero  ‘  order  values. 


Imuin  0  5 1'  1  ^  t  f  f  T* 
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4.3  Skin  Friction  and  Heat  Transfer 

To  boundary  layer  order,  the  ehear  Btreaa  at  the  surface  Is 


(85) 


Although  It  has  been  assumed  that  ^  order  to  provide  a 

better  fit  to  the  true  viscosity- temperature  behavior  over  the  entire  temperature 
range  encountered  In  a  hypersonic  boundary  layer,  this  approximation  Is  used 
only  In  the  differential  equations,  and  Is  not  used  to  evaluate  shear  and  heat 
transfer  at  the  wall. 

In  our  transformed  system,  thus. 

Me  approximate  by  where  we  have  assumed  that  a  linear 

fit  Is  possible  over  the  small  temperature  range  75-^  . 

/  t  £7^  (87)  whereTfe)  *  *7*  i  *^3 


y^~o 


The  density  at  y  ■  0  Is  given  by  the  perfect  gas  law, 

where  ^ 0  ) 


^(39)  -"/f  18 


th 


the  zero  order  density 
at  the  wall  . 


The  shear  stress  can  then  be  expanded  in  the  form 


If^  77  y/J 


11 


(90) 


since  jf  “-^Pr 


20 


The  skin  friction  coefficient,  defined  by 

is  given  by 

Cp^^(l*£Cf;)  (92) 


If  we  attempt  to  estimate  the  magnitude  of  this  slip  correction  to 


the  skin  friction  on  the  same  basis  that  we  estimated  the  correction  to  the 
induced  pressures,  we  can  approximate  the  resul  c  lor  ~t  /M  ifk  fu/  by 


(97) 

_  //^ 


(96) 
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For  this  case  of  the  adiabatic  wall,  it  is  possible  to  make  an  estimate  of  the 


magnitude  of  the  second  term  In  the  expression  for 


by  its  Blaslus  value 
o 


Then,  we  may  write 


11 


by  approximating 


(100) 


where 


L 


1 


(64) 


^  «  1. 


310  for  air 


If  we  approximate  Pj^  by  -1,  as  suggested  earlier,  the  result  is  then 

f,^^je  -.ietfi-l?/  (101) 


This  also  suggests  that  It  may  be  possible  to  predict  C_  from 

jf 

measured  values  of  which  are  obtained  from  the  measurements  of  induced 

pressures  on  adiabatic  plates.  If  the  above  procedure  represents  the  phenomenon 
adetjuately,  then  the  deviations  from  the  strong  interaction  pressures  could  be 
used  in  a  formula  of  the  form 

C I  -  I  (103) 

where  p^^  could  be  estimated  from  available  data. 
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The  heat  transfer  to  the  wall  may  be  estimated  more  readily,  but  it 
should  be  noted  that  the  energy  transferred  to  a  wall  In  slip  flow  Is  given  by 


(104) 


where  q  Is  the  Fourier  heat  conduction  term 
Is  the  slip  velocity 
Is  the  shear  at  the  wall. 

In  tezas  of  our  expansion  parameter,  £  , 

f.  -JiTi  - 


/?■  ftfo 


Since  W,  7*^  •  Eft'S  ft 


(106) 


The  total  energy  transfer  to  the  wall  Is  then 


(105) 


(10?) 


1 1 


For  any  wall  temperature  and  Pr  «  1, 

I.  ^  -.^2- 


0  * 


since 


3‘1‘  A  'J  'l^  ^ouyo. 

Thus,  for  Prandtl  No.  &  1,  the  slip  work  term  is  exactly  balanced  by  the  change 
in  the  Fourier  heat  conduction  to  the  wall,  and  there  is  no  net  change  in  the 
energy  transfer  to  the  wall. 
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5  SUIMMg  AND  CCaiCIJJSIOHE 


The  eolutloa  to  the  problem  of  the  effects  of  slip  boundary  conditions 
In  the  etrong  lateracvlon  region  hss  been  preaentsd,.  Since  the  appropriate 
expansion  parameter  la  of  ,  tl^  effects  of  slip,  as  obtained  from 

this  solution  are  only  one  of  several  possible  higher  order  effects. 

In  view  of  the  complexity  of  the  general  solutions,  crude  estimates 

of  the  slip  effects  on  Induced  pressures  and  skin  friction  were  made  only  for 

the  adiabatic  wall.  These  are  presented  In  Figures  1  and  2.  In  addition.  It 

was  found  that  for  the  case  of  arbitrary  constant  wall  temperature,  and  Fr  ■  1, 

fch 

there  la  no  change  in  the  heat  transfer  from  sero  order  strong  Interaction 
solutions* 

In  order  to  evaluate  the  correctness  of  the  estiostss  for  induced 
pressure  and  skin  friction,  It  would  be  necessary  to  integrate  several  functions 
of  the  sero^^  order  quantities*  It  would  be  desirable  to  do  so,  since  this 
would  show  whether  or  not  the  significant  departures  from  the  strong  interaction 
theory  in  the  induced  pressure  can  be  attributed  to  slip  at  the  surface*  If  it 
is  decided  that  slip  phanomana  are  present,  and  are  Important,  then  the  use  of 
the  ccuiplsts  solutier.2  would  aid  in  tb*  prediction  of  ekln  friction  from 
preeaure  data. 

Additional  work  Xi  obviously  necessary  on  the  general  problam  of  the 
hypersonic  leading  edge,  both  from  the  viewpoint  of  kinetic  theory,  starting 
from  the  tip,  and  also  from  the  viewpoint  of  e  set  of  continuum  conservation 
equations  such  as  the  Nsvler-Stokes  equations. 

The  apparent  Inadequacy  of  the  hypersonic  boundary  layer  equations 
In  predicting  aerodynamic  quantities  suggests  that  higher  order  effects  may  be 
both  present  and  Blgnlflcent, 
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ESTIMATED  INDUCED  PRESSURES  FOR  ADIABATIC  WALL 
WITH  FIRST  ORDER  SLIP  BOUNDARY  CONDITIONS 
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APPEMDIX  A 


EVALUATION  OF  SMALL  PARAMETER 


J.= 


(A-l) 

Vj.  Jy  =  (*-2) 

>  5^  3,  'A- 

»  /fc-  TSj  /t^V  7i^ 


(A-3) 


.-.  f '  (S  r5-f  t 

•  0  u  •  •  riio 


(A-4) 


(A-5) 


For  adiabatic  floM,  and  Fr 

£  •.XZ/HiSJN 

w  j . 


a-6) 
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APPENDIX  B 


The  expansion  of  the  purely  x-dependent  quantities 


'fy  r  [_!  ed')  -~J 

fPia 

(B-1) 

^  =  (f.  [  1  +£(S.  +  O^'-)  + 

(B-2) 

«het.-  /?= 

^  /f, 

W;-  UioLl+^*^'  1-0^')  ■>■■'] 

(B-3) 

f--  f,C\  "-J 

(B-4) 

vhere  the  coefficients  ,  '  »  ■  * 

are  constants 

To  first  order  In 

e 

/ 

from  (21) 

(B-5) 

S 

-3/,  k 

from  (33) 

(B-6) 

r» 

—  1 1 

M, 

from  * 

constant 

Therefore,  to  hypersonic  order 

t 

ft  %  K  (B-9) 

(A,  *  6*  .  (B-10) 
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APPENDIX  C 


The  evalxiatioii  of  the  arbitrary  constants  in  the  general  solution 
requires  two  additional  boundary  conditions,  since  the  general  solution 
identically  satisfies  the  conditions  that  as 

On  physical  grounds,  we  can  require  that  the  displacement  and  momentum  thick¬ 
nesses  be  finite,  since  otherwise  our  entire  approach,  being  based  strictly 
on  the  boundary  layer  hypothesis,  would  be  invalid.  This  requirement  is 
also  in  keeping  with  our  assumption  that  the  effects  of  first  order  slip  can 
be  expressed  as  a  perturbation  on  the  no  slip  solution.  In  addition,  this 
requirement  guarantees  that  the  only  singularity  in  the  vertical  velocity, 
and  hence  pressure,  will  occur  at  the  leading  edge. 

The  displacement  thickness,  in  hypersonic  flow,  can  be  written  as 

'’’H  T  »■> 

neglecting  terms  of 
-  f  ^  7 

Expanding  in  t.  yields  ^  Jo  -f  /  CC-i) 


where 


function  of  x 


(C-2) 


fi.inctlon  of  X 


(C-3) 
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<C-4) 


The  momantus!  thickness  Is 


e  '*■^1  ^ 


vhich  can  he  written  similarly  as 

©  =  6,*te, 


where 


function  of  x  (C-5) 

function  of  x  (C-6) 

» 

-f“  L  X  function  of  X  (C-7) 


The  requirement  of  finite  thickness  at  a  station  x  thus  requires 
that  the  inte«ala 

"*> ,  ^-Q.  >>“■" 


♦■V* 


all  be  bounded. 

Since  the  results  of  the  zero'""  order  solution  are  known  to  yield 
finite  values  for  (C-2)  and  (3-6),  we  are  then  led  to  the  condition  that 


and 

.tO 


be  finite. 


Thus  j  ^  ^  I 


<.  (JUf 


(C-8) 


(C-9) 
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since  the  eaymptotic  expansion  for  f  '  Is  (from  reference  [  19)) 


•H'  *«  - 

(C-IO) 


we  can  write 


^  ®  (C-ll) 

and  li  i  Si  \ 

<r  ezi 

the  nature  of  the  general  solutions  for  fj^’  reveals  that  the 
exponential  terns  in  the  integral  will  yield  only  a  finite  contrlbuttonc 
Therefore ,  we  may  write 


I  -2<') -^7! 


and 


/r^y 


<  ao 


(0-14) 


This  leads  to  the  requirement  that 


/  /  //  '^'5 1  ^ 

The  general  solution  for  is 


(C-16) 


7 


j  ^  — I  P  ,  where  F  Is  a  constant 


'  /4-w)  *  -^d’« 


i 


«/*)  (C-17) 
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or  since 


The  condition  for  finite 


as  -*®0j 

^  “s  exponential  terms  (C-18) 

7-r 

•or  finite  J  requires  that  B  •  0.  Thus 

«:-l9) 


(C-19) 


where  A  is  determined  by  the  temperature  jump  condition  at  k)  * 
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APPENDIX  C  “  (continued) 

The  requirement  that  J" fi  eLj  5  ^  be  finite  as  y  — enables 
us  to  determine  one  arbitrary  constant  in  the  solution  for  The  general 

solution  for  f^^'  is  in  the  form 

A“CC-^ 


by  the  behavior  as  . 

Knowledge  of  the  asymptotic  behavior  of  the  zero  order  solution  is 
essential  in  order  to  evaluate  0. 


Writing  ctj 

.  r  a  >-  a'  / 


e 


(C-20) 


since  the  Integral  ^ 


00 


is  finite. 


For  large  values  of 


yU. 


7 


[19  ] 


,  the  integral  can  be  expanded  asymptotically 

-■Z/F+' 

-  (C-21) 


I  0 


\  r  i  , 


/p-Oa J 


\v» 


I  ^  v- Zv,  r>j-rj)  X  ..Jc/ycc- 


24) 


(58) 
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(C-25) 


For  the  case  where  ^  ,  heat  transfer  cage 

*  ■■  ■ 

For  the  case  where  E  <■  0 ,  adiabatic  waL  1 

V*^  -ij  Voa  -  <Lj  *  •  0:-26> 

>1-40  I 

Since  the  escpainfllone  for  the  Integrals  are  onl^  c-f  exponential  order. 


and  contribute  finite  quantities, 


*  Voo  -  f  J 


(C-27) 


=  Von 


Similarly  for  the  expansion  for  larg«  can  be  ' 

V,  —  ■*■  '  -  <C  39) 


(C-28) 


be  written 


(C-30) 


Considering  the  leading  terms  in  the  expaaslon  for  f^  ^ 

u3 

Z'^'  I 

The  requirement  that  ^  ft  he  bounded  is  sa  tlsfled  only  be  setting 

the  arbitrary  constant  P  equal  to 

or  X)-  id  ,  (c-32) 


(C-31) 
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